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Abstract 

Reaction-diffusion systems which include processes of the form A+A — ► A 
or A + A — > are characterised by the appearance of 'imaginary' multiplica- 
tive noise terms in an effective Langevin-type description. However, if 'real' 
as well as 'imaginary' noise is present, then competition between the two 
could potentially lead to novel behaviour. We thus investigate the asymptotic 
properties of the following two 'mixed noise' reaction-diffusion systems. The 
first is a combination of the annihilation and scattering processes 2 A — > 0, 
2 A —> 2B, 2B — ► 2 A, and 2B — > 0. We demonstrate (to all orders in pertur- 
bation theory) that this system belongs to the same universality class as the 
single species annihilation reaction 2A — ► 0. Our second system consists of 
competing annihilation and fission processes, 2A — > and 2A — ► (n + 2)A, 
a model which exhibits a transition between active and absorbing phases. 



However, this transition and the active phase are not accessible to pertur- 
bative methods, as the field theory describing these reactions is shown to be 
non-renormalisable. This corresponds to the fact that there is no stationary 
state in the active phase, where the particle density diverges at finite times. 
We discuss the implications of our analysis for a recent study of another 
active / absorbing transition in a system with multiplicative noise. 

PACS Numbers: 02.50. -r, 05.40. +j, 82.20. -w. 
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1 Introduction 



Recently the effects of fluctuations in reaction-diffusion systems have attracted 
considerable attention (for reviews, see Refs. 0])- In sufficiently low spatial di- 
mensions the presence of microscopic particle density fluctuations causes traditional 
approaches, such as mean-field rate equations, to break down. This has led to the 
introduction of field-theoretic methods, based on 'Hamiltonian' representations of 
the associated classical master equation || |], || . These methods allow fluctuations 
to be handled in a systematic manner. The first system to be analysed in this way 
was the single species annihilation reaction A + A — > 0, || 0] where it was shown 
with renormalisation-group (RG) methods that for dimensions d < 2, the average 
density n(t) decays to zero at large times according to the power law 

W)~E d r d '\ (l) 

with Ed denoting a universal amplitude (for uncorrelated initial conditions), while 
n(t) ~ E 2 t~ 1 lnt in d = 2. Furthermore Peliti has demonstrated that the coag- 
ulation reaction A + A — > A belongs to the same universality class as the pure 
annihilation process A + A — > |J. Physically the anomalously slow decay of 
eq. (HD results from the emtzcorrelation of particles in low dimensions. Due to 
the 'reentrancy' property of random walks for d < 2, once two particles are in 
close proximity they will then tend to react rather quickly. Hence at large times 
the remaining (unreacted) particles are likely to be situated far from their nearest 
neighbours (i.e. the particles become anticorrelated). 

Markedly more complex behaviour may arise once particle production processes 
are also permitted. For example, the 'Branching and Annihilating Random Walk' 
(BARW) system defined by the reactions 2A — > and A — ► (m + 1)A, displays a 
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dynamic phase transition between an 'active' (n(t) — > n s > for t — > oo) and an 



'inactive and absorbing' state (n(t) — ► for t — ► oo), with a remarkable difference 
between the cases of odd and even number of offspring m |,|. For odd m and 
d < 2 the transition is basically characterised by the critical exponents of directed 
percolation (DP) [|K], [11], |12|, 0, whereas for even m and d < d' c ^ 4/3 the phase 
transition is described by a new universality class, with the density in the entire 
absorbing phase decaying according to the power law in eq. (P |8|, |9|, [14 . 



A powerful method for the analysis of such systems is provided by the RG 
improved perturbation expansion |7|, [TJ|, [ljj . However, once a field-theoretic action 
for the system has been derived (from a microscopic master equation), it is also 
possible to write down effective Langevin-type equations, where the form of the 
noise can now be specified precisely, without any recourse to assumptions and 
approximations f2j. The nature of the noise can look somewhat peculiar in this 
representation, for example in the A + A — > reaction we have the exact equation 
for the field a(x, t): 

d t a(x, t) = D V 2 a(x, t) - 2A a(x, t) 2 + a(x, t) rj(x, t) , (2) 

where D is the diffusion constant, A the reaction rate, and 

( v (x,t)) = 0, (7](x,t)r](x',t')} = -2X5 d (x - x')5(t - t') . (3) 

Hence the noise rj is imaginary, a rather counterintuitive result. Recently Grin- 
stein, Munoz and Tu [|17]] have studied an equation superficially similar to that 
given above, with the aim to model active / absorbing transitions in autocatalytic 
chemical processes. In the special case of a scalar field with a quadratic nonlinear- 
ity, their model is defined by the equation 

d t a(x, t) = D V 2 a(x, t) — r a(x, t) — u a(x, t) 2 + a(x, t) r](x, t) , (4) 



where 

<77(a:,t)> = , (ri(x,t)ri(x',t')) = 2u5 d (x -x')6(t-t') . (5) 
Notice that the noise in eq. © has the opposite sign to that considered previously 



(i.e. the noise of Ref. |T7j is real). However, this is an important point since 
a positive sign in the noise correlator leads to divergences in the renormalised 
parameters of the theory: 

v R =Zv i f d d k 1 

with Z = — and IAr) = / -, — r-r ttttt • (6) 

„ \-vUr) y ' J (2n) d r + Dk 2 y 1 

Ur = Zu v ; v ; 

Hence new singularities emerge when the denominator of Z vanishes. These diver- 
gences have not been present in earlier field-theoretic studies of reaction-diffusion 



systems. Certainly if the noise in the model of Ref. [17fl were resulting from the 
reaction A + A — ► 0, its correlator should have a negative sign, as described above 
(and hence no extra divergences would appear, rendering much of the interest- 
ing behaviour in Ref. |T7| obsolete). It is therefore not clear to us how the real 
multiplicative noise of eqs. (|),(§P could be the only type generated — we believe 
that internal, imaginary reaction noise should generically be present as well. Con- 
sequently the physical mechanism behind the noise analysed in Ref. [H] remains 
somewhat obscure. 

One of the objects of this letter is to see if equations similar to that anal- 
ysed by Grinstein, Munoz and Tu (with real noise, and hence potentially novel 
behaviour) can be derived consistently for certain reaction-diffusion systems using 
field-theoretic methods. Potentially at least, in an emerging competition between 
'real' and 'imaginary' noise contributions, the 'real' component might prevail in 
certain circumstances, conceivably leading to the scenario discussed in Ref. [ [T7|| . 
However, our main finding here is that although we have analysed two systems 



where the noise correlator has both positive and negative components, we have 



been unable to recover the new features discussed in Ref. ||17|| . In fact, in the first 
of our model systems, a combination of two-particle annihilation and scattering 
processes for the species A and B, the ('imaginary') reaction noise dominates the 
long-time behaviour, which is described by the asymptotic power law ([]]). Neither 
can the 'imaginary' noise terms be neglected in our second reaction-diffusion sys- 
tem, namely combined annihilation and fission processes of a single species A. In 
this system a perturbative RG analysis breaks down in the active phase and at 
the dynamical phase transition separating it from the inactive state. Therefore al- 
though we can only address the inactive phase, which is again governed by the pure 
annihilation model, we believe this system cannot reproduce any of the features in 
Ref. HTfl. 



In the following Sec. 2, we present a more thorough discussion of how 'imaginary' 
noise terms emerge in processes dominated by two-particle annihilation reactions. 
On the other hand, problems belonging to the directed-percolation (DP) univer- 
sality class, as described by Reggeon field theory |H| [11], [12], [13], can be faithfully 
represented by a simple Langevin equation for the local particle density with 'real' 
noise. In Sec. 3, we next present our scattering / annihilation model, which we first 
analyse to one-loop order, and then to all orders in perturbation theory by solving 
the coupled Bethe-Salpeter equations for the vertices. In Sec. 4, we proceed to dis- 
cuss the annihilation / fission reaction system, and show that while its properties 
in the inactive state may be analysed using field-theoretic methods, this is not the 
case in the active phase or at the dynamic transition itself. Finally, we summarise 
and discuss our results in the light of the recently proposed transition scenario for 



'real' multiplicative noise problems fu |. 
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2 'Real' vs 'imaginary' noise in reaction-diffusion 

systems 

Before turning to our investigation of models with competing 'real' and 'imaginary' 
noise terms, we briefly outline and review the general issue of how to systematically 
include fluctuation effects in reaction-diffusion systems. Above the upper critical 
dimension, a qualitatively correct analysis may be obtained from the associated 
mean-field rate equations for the average particle densities. Below this dimension, 
where fluctuations become important, it is tempting to apply a Langevin equa- 
tion approach, motivated by the success of this technique in equilibrium critical 
dynamics. However, one has to be aware that these typically irreversible reactions 
constitute a dynamical system far away from thermal equilibrium. Thus there is 
no fluctuation-dissipation theorem available which could serve as a guide to the 
appropriate form of the Langevin noise correlations. One could of course just try 
the simplest ansatz, namely some form of white noise multiplicatively coupled to a 
certain power of the particle densities, in order to ensure that all fluctuations van- 
ish when there are no particles left (i.e., in the absorbing state). But, as we shall 
see shortly, at least for processes dominated by two-particle annihilation reactions, 
this generically leads to an incorrect analysis. 

Thus, in order to systematically include the effects of microscopic density fluc- 
tuations in low dimensions, one can instead start with the corresponding classical 
master equation, then represent this stochastic process by the action of second- 
quantised bosonic operators, and finally use a coherent-state path integral repre- 
sentation to map this system onto a field theory. This mapping itself is a standard 
procedure, and is described in detail in Refs. @, |], [5], 0, for instance. Apart from 
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the continuum limit that is usually taken, this procedure provides an exact mapping 
of the initial master equation, and involves no assumptions whatsoever regarding 
the form of the noise, the relevance or irrelevance of certain terms etc. Note that 
the resulting bosonic theory applies only to systems where there is no restriction 
on the particle occupation number in the microscopic model. For the description 
of exclusion processes where the site occupation numbers are restricted to or 1, 
obviously a fermionic representation is more useful. 

For example, for the simple two-particle annihilation reaction A + A —>■ the 
ensuing field-theoretic action reads 

S = Jd d x Jdt [a(d t - DV 2 )a - A(l - d 2 )a 2 ] , (7) 

where we omit boundary terms relating to the initial conditions and the projec- 
tion state. Here D denotes the diffusion constant, A the annihilation rate, and 
a(x,t) and a(x,t) are bosonic fields. The stationarity conditions ('classical field 
equations') 5S/5a = and 5S/5a = yield, respectively, a = 1 and the mean-field 
rate equation dta = DV 2 a — 2Aa 2 . Thus, within the mean-field approximation 
we can identify a(x,t) with the local 'coarse-grained' particle density. However, 
fluctuations in a(x, t) may not be simply related to density variations, as can be 
seen by performing the shift a = 1 + a, 

S = Jd d x Jdt [a(d t - DV 2 )a + 2\da 2 + AaV] . (8) 

Integrating out the 'response' field a in the functional integral JVaVdexp(-S) 
then leads precisely to eq. (0) with the negative noise correlator (|3]). Physically 
this counterintuitive result corresponds to the antzcorrelation of particles in low 
dimensions. Furthermore, power counting reveals that this noise, which originates 
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from the quartic term in the above action, becomes relevant below a critical dimen- 
sion d c = 2. Because of this pure imaginary noise, a(x,t) clearly cannot represent 



a physical density variable. Although {a(x, t)) is equal to the mean density n(x : t), 
similar relations do not hold for higher correlators, see Ref. ||. Moreover, since 
a(x, t) is not the density field, this means that the noise must also take a non- 
standard from — in fact rj in eq. @ represents only the contribution to the overall 
noise from the reaction process. In reality, of course, diffusive and reaction noise 
can never be disentangled from one another, and consequently there is no partic- 
ular reason why the noise in eq. (|3|) should be 'real', i.e., described by a strictly 
positive correlator. 

However, it is possible to give descriptions where the diffusive noise does appear 
explicitly. One approach begins with a Langevin equation including both real 
reaction noise and diffusive noise. In that case the field a(x,t) represents a coarse- 
grained local density. This is the approach taken, for example, by Janssen in 
Ref. fl2fl . A second possibility is to begin with the representation (0) and then 
obtain an equivalent description in terms of 'density' variables by considering the 
canonical transformation a = exp(— p)p, a = exp(p), da = p (see Ref. |T3fl). The 
resulting effective action in terms of the new 'density' fields p and p then includes 
a term oc —p 2 p 2 which corresponds to pure real noise. However, the ensuing field 
theory contains an extra 'diffusion noise' contribution. 

Hence we see that the 'naive' Langevin equation (i.e., eq. (|2]), but with a positive 
noise correlator and no diffusion noise) does not provide an appropriate effective 
description of the above system. On the other hand, for the standard Gribov 



process A — > 0, A — > A + A, A + A — > A [|T|, |13fl , the action in terms of the shifted 



9 



fields reads 

S = Jd d x Jdt \a[d t + D(r - V 2 )]a - aa 2 a + Xaa 2 + AaV] . (9) 

Here D and A represent the diffusion constant and annihilation rate as before, a 
is the branching rate, and r = (fi — o)/D where \i denotes the spontaneous decay 
rate. Integrating out the response fields now yields 

dta(x , t) = D (V — r)a(x, t) - \a(x,t) 2 + r](x,t) , (10) 

with 

(i](x, t)) = , (i](x, t)r](x', t')) = 2[a a(x, t)- \a(x,t) 2 }6 d (x-x')6(t-t') . (11) 

It should be noted, however, that the effective coupling entering the perturbation 
expansion is actually oc crA, which leads to an upper critical dimension d c = 4. 
Consequently, the term Xd 2 a 2 becomes irrelevant, and we are left with a pure 
positive definite noise correlator oc aa(x,t) (or 'square-root' multiplicative noise, 
if we replace rj with y/arf). After a simple rescaling, the action is readily mapped 
onto Reggeon field theory for directed percolation |10, [11, [1^, |j~3| . Thus, here we 



encounter the generic case where a(x,t) can be identified with a coarse-grained 
density field. Physically, the above reactions lead to particle clustering, and local 
densities indeed constitute a natural choice for the order parameter field. We also 
remark that the signs and magnitudes of the prefactors (which may even be chosen 
to be imaginary) of the cubic nonlinearities in the action (§) do not matter provided 
their product — o~\ remains real and negative. 

In the above analysis, we have brought out the contrasting features of two cate- 
gories of reaction-diffusion systems — containing either imaginary (anticorrelating) 
or real (clustering) noise. However, reaction-diffusion systems with a multiplicative 
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noise term oc arj are characteristically governed by pair reaction processes and thus 



have 'imaginary' noise. This is in contrast to the ansatz in Ref. fl7 |. Of course, 
one might argue that if there are both 'real' and 'imaginary' noise contributions 
present, then the 'real' parts might prevail and lead to the scenario discussed in 



17]. This possibility motivates the following two case studies of combined scatter- 



ing / annihilation and annihilation / fission reactions to which we now turn our 
attention. 



3 The scattering and annihilation process 

The first react ion- diffusion system we want to consider consists of the four reaction 
processes 

A + A^(&, A + B + B , B + B -> A + A , and B + B -> , (12) 

which occur at rates Xaa, ^ab, ^ba, and Xbb, respectively, and with diffusion 
constants Da and Db for the A and B particles. We choose uncorrelated initial 
conditions where the A and B particles are distributed randomly. Physically the 
above reaction scheme might occur if the A particles could undergo a scattering 
process turning into B particles, and vice versa, in addition to the presence of 
the annihilation reactions. In order to systematically include the effects of mi- 
croscopic density fluctuations in low dimensions, we represent the corresponding 
master equation by a coherent-state path integral (see Sec. 2). In terms of the 
continuous fields a, a, b, b, the diffusivities Di ^ 0, continuum reaction rates {A^}, 
and the initial homogeneous densities rii (where i,j = A,B), the action reads (for 
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S = Jd d x Jdt [a(d t - D A V 2 )a + b{d t - D B V 2 )b + 2X AA aa 2 + X AA a 2 a 2 + 
+2X BB bb 2 + X BB b 2 b 2 + 2X AB da 2 + X AB d 2 a 2 — 2X AB ba 2 — X AB b 2 a 2 + 

+2X BA bb 2 + X BA b 2 b 2 -2X BA ab 2 - X BA a 2 b 2 -n A a5(t) -n B b5(t)] . (13) 

If we now integrate out the response fields a and b from the functional integral 
J T>aDdT>b'Dbexp(—S), we find that the above reaction-diffusion system can be 
described exactly by a pair of Langevin-type equations 

d t a(x, t) = D A V 2 a(x, t) - 2(X AA + X AB ) a(x, tf + 2X BA b(x, tf + rj A {x, t) , 
d t b(x, t) = D B V 2 b(x, t) - 2(X BB + X BA ) b{x, t) 2 + 2X AB a{x, t) 2 + r] B (x, t) ,(14) 

with noise correlations 

( VA {x,t)) = ( VB {x,t)) = , (15) 
(Va{x, t)r] A (x', t')) = [X BA b{x, t) 2 - (X AA + X AB ) a{x, t) 2 ] 5 d {x - x')8{t - t') , 
(Vb(x, t) VB (x', t')) = [X AB a(x, t) 2 - (X BB + X BA ) b(x, t) 2 ] 5 d (x - x')5(t - t') . 

Hence, as desired, we have constructed a system where in a Langevin-type for- 
malism we have terms of both signs present in the correlator. Thus we can now 
attempt to answer the question of whether this 'competition' alters the structure 
of the theory in low dimensions where fluctuations are of vital importance. 

Power counting on the action fll~3|) reveals that all the reaction rates {A^} have 
dimension ~ Ai 2_a! , where \x denotes a momentum scale. Hence we expect to find 
a critical dimension d c = 2, below which fluctuations change the mean-field be- 
haviour qualitatively and the theory must be renormalised. As in the pure annihi- 
lation model (||), this renormalisation is simple since the diagrammatic structure of 
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Figure 1: The temporally extended vertex function Aaa(^, s) to one-loop order. 

the theory does not permit any dressing of the propagators. Hence the only renor- 
malisation required is that for the reaction rates. We now define \aa = ^aa + ^ab, 
^ab = Aab, ^ba = ^ba and \ BB = \ BB + \ba- The temporally extended vertex 
function for \ A A{k, s) to one-loop order is given by the sum of diagrams shown in 
Fig. 1 (here s is the Laplace transformed time variable; time runs from right to 
left). The diagrams for the other vertex functions look quite similar. Evaluation of 
these one-loop diagrams yields the following form of the renormalised reaction rates: 
9ij = CA^k, s)| fe2/4=M2iS=0 /A/A where e = 2 - d, and C d = T(2 - d/2)/2 d - 1 n d / 2 
is a geometric factor. This leads in a straightforward manner to the following 
one-loop RG beta functions = /idg^/d/i: 

Paa = 9aa{-£ + 9aa) + 9ab 9b a , (16) 

Pab = 9ab{-£ + 9aa + 9bb) , (17) 

Pba = 9ba{-£ + 9aa + 9bb) , (18) 

Pbb = 9bb(-^ + 9bb) + 9ab 9ba ■ (19) 

In fact it can be shown that the above one-loop beta functions are actually 
exact to all orders in perturbation theory. This is readily accomplished by writing 
down the full coupled Bethe-Salpeter equations for the vertices. Diagrammatically, 
this corresponds to replacing either the right-hand or the left-hand bare vertices in 
all the one-loop contributions (see Fig. 1) by their fully renormalised counterparts. 
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This freedom of choice immediately implies the relation 

X AB (k, s)/X AB = X BA (k, s)/X BA = N(k, s) . (20) 

After absorbing the diffusivities Di into the bare couplings A^- and the full ver- 
tex functions Xij(k, s), respectively, and introducing the abbreviation Ii(k,s) = 
(2ir)~ d J d d p [p 2 + [k 2 /A) + (s/2D i )]~ 1 , the coupled exact Bethe-Salpeter equations 
can be explicitly written as 

X AA (k, s) [1 + X AA I A (k, s)] + X AB (k, s) X BA I B (k, s) = \ AA , (21) 
X AA {k, s) X AB I A (k, s) + X AB (k, s) [1 + X BB I B (k, s)) = X AB , (22) 

together with a second pair of equations which follow by interchanging A <-> B in 
eqs. (|2l| ) and (^2|). These coupled linear equations (|20|) — C|22[) for Xj,j(k, s) are solved 
by 

N(k, s)- 1 = [1 + X AA I A (k, s)] [1 + X BB I B (k, s)} - X AB X BA I A (k, s) I B (k, s) , (23) 
and 

[X AA (k, s)/X AA N(k, s)] = l+ [X BB {1 - X AB X BA /X AA X BB )} I B {k, s) , (24) 
[X BB {k, s)/X BB N{k, s)] = l + [X AA (1 - X AB X BA /~X AA X BB )} I A (k, s) . (25) 

At the normalisation point one has Ij(2/i, 0) = CdH~ e /e, and after some tedious but 
straightforward algebra eqs. (|23| ) - (|25|) yield again the beta functions fll6l)-(|l9D. 



We can now examine the above eqs. (pT6|)— ( IE ) , which we have just demonstrated 
to hold to all orders in perturbation theory, for fixed point solutions defined by 
f3ij{{g*j}) = 0. For d > 2 we find, as expected, merely the trivial Gaussian fixed 



point where all g*j = 0. However, for d < 2, the only stable fixed points are those 
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describing uncoupled annihilation processes, i.e., 

9*AB = 9*BA = , g*AA = 9*BB = e ■ ( 2 6) 

Furthermore, there are also other solutions, for example the fixed line 
< 

c = 9*ab9*ba — e2 /4 fixed but arbitrary, and 



2g AA = e±Ve^4c, 2g* BB = e T Ve^4c ; (27) 

but these, like the Gaussian fixed point, turn out to be unstable for d < 2. 

Hence the above annihilation / scattering model (|T2| ) asymptotically becomes 
rather simple, and in fact lies in the same universality class as single-species an- 
nihilation (with respect to both the decay exponent and amplitude). Hence each 
species of particle decays according to eq. (P as t — > oo for d < 2. Physically this is 
a result of the 'reentrancy' property of random walks — as soon as two particles are 
in close proximity, they will rapidly annihilate, even in the presence of scattering 
processes. Therefore we conclude that, for this system, the presence of 'real' as well 
as 'imaginary' noise has not introduced any novel behaviour. We finally remark 
that the above results also apply in the extreme asymmetric situation where, say, 
^ba — but Xab > originally, i.e., when there is spontaneous unidirectional 
transformation of pairs of A particles into pairs of B particles, but not vice versa. 
At least in this special case, our result that this pairwise transmutation is irrele- 
vant in the long-time limit, appears nontrivial. For example, in the related case of 
DP processes for A, B particles with coinciding critical points, which are coupled 
via the reaction A — > B, the usual DP critical exponent (3 is replaced by a much 



smaller density exponent as a consequence of an ensuing multicritical point [|T^] . 
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4 The annihilation and fission process 



Our second react ion- diffusion system consists of the processes 



A + A 



and A + A -> (n + 2)A 



(28) 



to which we assign the annihilation rate A and 'fission' rate a n . Note that these 
processes differ from the 'Branching and Annihilating Random Walks' || |9], |14f 
mentioned earlier in that offspring particles can only be produced upon collision 
of two A particles. The corresponding action derived from the master equation 
describing the reactions (|28|) reads in terms of the unshifted continuous fields a(x, t) 
and a(x, t) 



S= d d x dt 



d(d t - DV 2 )a - A(l - a 2 )a 2 + a n (l - a n )d 2 a 2 



(29) 



(the terms depending on the homogeneous, uncorrected initial density distribution 
and on the projection state have been omitted here). Once again we point out 
that this theory is valid only for unrestricted particle occupation numbers in the 
microscopic model. It is quite possible that altering the microscopic rules for site 
occupancy (for example by allowing only or 1 particles at a site) may change some 



of our later conclusions |19| . If we now proceed by performing the shift a — 1 + a, 
the effective action becomes 



S 



d d x dt 



a(dt - DV 2 )a + (2A - na n ) da 2 + 



+ A — 



n{n + 3) 



n+2 



1=3 



n + 



a l a 2 



(30) 



2 ~ ' " ^ v / 

//all vertices a 1 a 2 for / > 3 are neglected, this field theory becomes equivalent to a 
nonlinear Langevin equation 



d t a(x, t) = D V 2 a(x, t) + [na n — 2A) a(x, t) 2 + a(x, t) n(x, t) 



(31) 
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(r](x,t)) = 0, (rj(x,t)r](x',t')) = [n(n + 3)a n - 2A] 5 d (x - x')5(t - t') , (32) 

which again describes competition between 'real' noise (associated with <j n ) and 
'imaginary' noise (associated with A). Upon comparing with the model of Ref. [ |17] , 



eqs. (£|), (|5|), we see that the annihilation / fission process apparently corresponds 
to their parameters r = 0, u = 2A — na n , and v = n{n + 3)<T n /2 — A. 

Notice that it is potentially dangerous to perform the shift a = 1 + a and 
then to arbitrarily omit certain nonlinearities || 14], due to the discrete symmetry 



a — > —a, a — *> —a, under which the action (|29| ) is invariant (for n even). This 
symmetry corresponds to local particle number conservation modulo 2, which is 
lost in the Langevin description based on (|30D . Furthermore, the neglected terms 
have the same scaling dimension as those retained. We therefore proceed with the 
analysis of the unshifted theory (^). We find that both the annihilation and fission 
rate have identical scaling dimension ~ /i 2_d and thus the upper critical dimension 
is again expected to be d c = 2. For d > 2, a description given by the mean-field 



equation (i.e. eq. ([H]) without noise) should become qualitatively correct. For 
na n < 2 A this leads asymptotically to a density decay n{t) ~ t^ 1 (with reduced 
annihilation rate \r = A — na n /2); for na n > 2 A, on the other hand, the density 
grows rapidly and diverges at t c = l/(na — 2X)n , where n is the initial density. 
Thus, there is no stationary state in the active phase. 

We now consider a one- loop analysis of the action ( p9[ ) for d < 2. To this order 
all the couplings associated with the interaction vertices in (^) are renormalised. 
For example the coupling a n is renormalised by the diagrams shown in Fig. 2a. 
However, in addition to this, other processes are also generated at this order: 2A — > 
nA (by a combination of fission and annihilation, see Fig. 2b) and 2A — > 2(n+ 1)A 
(by two successive fission reactions, see Fig. 2c). Furthermore, the mechanisms 
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Figure 2: One- loop diagrams for (a) the renormalisation of the a n coupling, (b) 
the generation of the process 2A — > nA from a combination of 2A — > (n + 2) A and 
2A — > 0, and (c) the generation of the process 2A — > 2(n + 1)A by two successive 
fission reactions. 

producing these processes (with increasingly large n) are not simply restricted to 
the one-loop level - - new particle creation vertices are effectively generated at 
each successive order in perturbation theory. Therefore the number of (relevant) 
higher-order couplings, each of which requires renormalisation, increases without 
bound as higher and higher orders of perturbation theory are considered. Hence 
we must conclude that the field theory is non-renormalisable: an infinite number 
of renormalisations would be needed to render the theory free from divergences. 

Nevertheless some further progress is possible by considering the original master 
equation for the annihilation-fission process ||20|| . Omitting the diffusive terms, we 
have 

dP ^ }]t) = A lat ^[(m J+ 2)K + l)P({. . . , m i+ 2, . . .}; t)-m J (m,-l)P({m J }; t)} 
+ a^Y,[(mi-n){m l -n-l)P{{. . . . . .};t)-m i (m i -l)P({m i };t)], (33) 

i 
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where P({mj};t) is the configuration probability for finding occupation numbers 
{mj at time t, and where A lat and crjf* are the lattice annihilation and fission rates, 
respectively. Using the relation m{t) = J2{ mi } TniP({mi}; t), eq. (|33| ) implies that 

^& = (na 1 :* - 2A lat ) m(m - 1). (34) 

Since rrii = 0,1,2,.. ., we see that m{m — 1) is non-negative, and hence that nerjf* = 
2A lat marks the transition point between the active and inactive phases. Note that 
actually at the transition {na 1 ^ = 2A lat ), the average density will remain constant, 
whereas in the active phase it will diverge. These conclusions can be confirmed 
by studying the shifted action (p0|). If we have the equality 2 A = na n for the bare 
field-theoretic parameters, then the bare cubic coupling vanishes. However, the 
structure of the higher-order vertices ensures that a cubic coupling cannot then 
be regenerated at any order in perturbation theory. Hence we can again conclude 
that 2A = na n is the transition point between the active and inactive phases. 
Furthermore in the inactive phase, where the annihilation mechanism dominates, 
and the successive generation of an infinite series of fission processes is probably 
suppressed, we might expect the density to decay as t~ d l 2 (for d < 2) due to the 
strong particle anfeorrelations which emerge as a result of the annihilation process 
in low dimensions. 

However, the non-renormalisability of the field theory means that we are unable 
to fully address the properties of either the active phase or the active / inactive 
transition. This failure may be associated with the fact that the active phase is 
not in a stationary state (at least in mean-field theory, and with unrestricted site 
occupancy m« = 0,1,..., oo, the density in this phase diverges in finite time). 
This behaviour, taken together with the massless nature of our field theory (i.e., 
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there is no term proportional to the field a(x,t) in eq. (|5ID), implies that the 
transition from the absorbing to the active phase cannot be in the DP universality 
class. Rather the transition is closer to being 'first order', as suggested by an exact 
evaluation of the correlation function at criticality |20]. As field-theoretic methods 
clearly cannot shed any further light on this problem, we hope that our analysis 
will stimulate further work using, for example, exact one-dimensional methods. In 
addition, numerical simulations presently in progress seem to indicate that this 
annihilation / fission system may display remarkably rich behaviour |I9|. 



5 Summary 

In this paper we have studied the effects of various types of noise in diffusion 
limited reactions. In section 2 we emphasised that 'naive' Langevin equations 
(with positive noise correlators) fail to accurately describe systems controlled by 
pair reaction processes, where the noise is in fact 'imaginary'. Physically this failure 
is associated with the anticorrelation of particles in low dimensions. On the other 
hand such a naive approach does indeed work for the Gribov process, where the 
noise turns out to be 'real' (related to particle clustering). 

We then studied two diffusion-limited reaction systems with both real and imag- 
inary noise components: the annihilation / scattering processes (O) and the an- 
nihilation / fission processes (j28D . We have shown (to all orders in perturbation 
theory) that the first of these belongs to the same universality class as the pure 
annihilation model in dimensions d < 2, while for d > 2 the mean- field rate equa- 
tions apply. However, the second system displays a transition between an active 
and absorbing state, which is not accessible to perturbative analysis. In both cases, 
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despite the competition between 'real' and 'imaginary' noise, we have been unable 



to recover any of the interesting behaviour discussed in Ref. |[L7| . In fact, consider- 
ing that the processes discussed here, along with BARW, are amongst the simplest 
reactions leading to both 'real' and 'imaginary' multiplicative noise, it is rather 
unclear which physical system might be described by the Langevin equation (|j) 
with purely 'real' multiplicative noise (||]), and thus display the nontrivial effects of 
Ref. (17|. 
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